Unbiased bases (Hadamards) for 6-level systems: Four ways from Fourier 
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In quantum mechanics some properties are maximally incompatible, such as the position and 
momentum of a particle or the vertical and horizontal projections of a 2-level spin. Given any 
definite state of one property the other property is completely random, or unbiased. For A'^-level 
systems, the 6-level ones are the smallest for which a tomographically efficient set of A'^ + 1 mutually 
unbiased bases (MUBs) has not been found. To facilitate the search, we numerically extend the 
classification of unbiased bases, or Hadamards, by incrementally adjusting relative phases in a 
standard basis. We consider the non-unitarity caused by small adjustments with a second order 
Taylor expansion, and choose incremental steps within the 4-dimensional nuUspace of the curvature. 
In this way we prescribe a numerical integration of a 4-parameter set of Hadamards of order 6. 

PACS numbers: 02.10.Yn, 02.60.Cb, 02.60.Pn, 03.65.Ta, 03.65.Wj, 03.67.Dd 



I. INTRODUCTION 

Mutually unbiased bases (MUBs) play an impor- 
tant role in quantum physics. Typically they are con- 
structed from the eigenstates of maximally incompatible 
properties.^ For example, in cases ranging from quantum 
optics to electronics, the state p of a photon polariza- 
tion or electron spin is often defined by a Stokes or Bloch 
vector ((?) whose components are the expectation values 
of the Pauli operators: p = {I + {a) ■ (t)/2. The Pauh 
operators a^, <Jy, Oz are incompatible because, given an 
eigenstate of any one of them, the outcomes of the oth- 
ers are completely random. Their eigenbases are said to 
be mutually unbiased. In general, two bases \vm) and 
\wri) are unbiased when all probabilities between them 
are equal: 



K^mluin)!^ = for all m,n 



(1) 



with A'^ as a normalization constant. 

MUBs are also manifest in Heisenberg's Uncertainty 
Principle: a definite momentum \p) implies a random 
position |a;), corresponding to its wavefunction of cqui- 
modular amplitudes 
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jpx/n 



(2) 



For a discrete system having N orthogonal states, or 
levels, a set of + 1 MUBs gives a sufRcienli^ and 
optimal^ way to determine its possibly-mixed state p. 
MUBs are thus useful for quantum state tomography. 
They also play a role in quantum cryptography,"* oper- 
ator algebraSf^ and Lie algebras<^ For a 2-level system, 
or qubit, the Pauli operators provide a useful set of 3 
MUBs. More generally, sets oi N + 1 MUBS have been 
constructed for systems whose number of levels A" is a 
prime number— or a power of a prime numberi^ In any 
case it is never possible to have more than A^-|- 1 MUBsii 

For a 6-level system, such as a spin-5/2 or a qubit- 
qutrit pair, there is a widespread and ongoing search for 
sets of MUBs4 6 is the smallest integer that is neither a 



prime nor a power of a prime, and, despite the qualitative 
and physical appeal of sets of A^-|- 1 MUBS, no more than 
3 MUBs have been found^ d°d^ for the 6-level systems. 

Most searches begin with the classification of unbiased 
bases represented by HadamardsiS. U that are unitary 
and comprise, in a standard basis |si), equi- modular ele- 
ments: 



U^U = I and 



(3) 



The columns of {si\U\sra) form an orthonormal basis of 
vectors \um) = U \sm) unbiased with respect to the stan- 
dard basis; we are using the term "unbiased basis" rel- 
ative to a standard basis. Quantum mechanics allows a 
choice of overall phase for each \si) and \um) which we 
use to fix, in the first row and column, (j)i„i — 4>ii — 0. 
Any reordering of the standard and unbiased bases' vec- 
tors (the permuting of rows and columns in U) is also 
immaterial. 

The classification of these equivalent Hadamards is an 
open question for A > 5^^ For A^ = 6 the known fami- 
lies of Hadamards have up to 2 parameters and are topo- 
logically connectedji^ it has been conjectured there are 
4-dimensional familie^ which would vastly expand the 
searchable candidates for sets of MUBs. The conjecture 
stems from an upper bound, know as the defect, on the 
dimensionality of any analytic set of Hadamards.-'^ 

Here we support the conjecture by numerical integra- 
tion of Hadamards of order 6 in 4 directions within the 
25-dimensional space of relative phases (j)im {} ^ ^ ^ 
to). We do this by Taylor-expanding to second order 
the non-unitarity of 6 x 6 matrices of equi-modular el- 
ements around known- and found-Hadamards. Neigh- 
boring Hadamards lie nearby in the paraboloid valley of 
the expansion. Apart from Tao's matrix,— we always 
find the valley, i.e. the nuUspacc of the curvature, to be 
4-dimensional, consistent with the defect, although that 
nuUspace can vary from Hadamard to Hadamard. By 
taking small steps constrained to the changing nuUspace, 
we numerically integrate Hadamards according to 4 dis- 
tinct parameters. 



2 



II. FOURIER AND OTHER KNOWN FAMILIES 



A common Hadamard is the Fourier matrix 



(sA F \Sm) = 



^2^tU-l){rn-l)/N 



N 



(4) 



from which there originate two 2-parameter afSne fami- 
hes. The first Fourier family, Fq, has components 



/II 111 1 \ 

1 qzi q'^Z2 q'^Zi q^Z2 

1 q^ q^ I q^ q^ 

1 q^zi Z2 q' 
1 q4 q2 ^ 

\ 1 q^zi q^Z2 q^ q^Zi qz2 J 



zi 



q^z2 

q^ 



(5) 



with q = e^^^*/^ and with Zi = e'"^! and Z2 = e*'^^ set- 
ting certain phases according to the affine parameters 
and <j)2- It is Hadamard because it comprises equi- 
modular elements and is always unitary — its columns 
are an unbiased basis (with respect to the standard ba- 
sis). Its transpose Fj {(j)^^ (^4), with (^3 and 4)i taking the 
(transposed) places of and (^2, is the second Fourier 
family. 

There are other Hadamard families: Beauchamp's and 
Nicoara's B^ijj)^^ which interpolates from Bjorck's cir- 
culant matrix Cg^^ and its conjugate to Diva's Dq{6)-}^' 
and Matolcsi's and SzoUosi's Mq{x) which connects Dq 
with Fq and F^ There is also Tao's matrix 5*6 com- 
prising 3rd roots of unityji^ which is disconnected from 
the resti^ 



III. NON-UNITARITY TO SECOND ORDER 



excluding the (pim = 4>ii = in the matrices' first row 
and column. 

The difficulty is in choosing the phases to ensure a ma- 
trix is unitary. We therefore define a measure of its non- 
unitarity as proportional to the sum of the probabilities 
between non-trivial pairings of its column vectors, 



/=(7V/2) l(e™|e„)l' >0, 



(7) 



which vanishes only when E is unitary (and thus a 
Hadamard). Taking the inner product (em|e„) and its 
conjugate in the standard basis, we find the non-unitarity 
/ is a simple function of the relative phases: 



/=EE 



cos 



n>m j—i 



Our methods will rely on Taylor-expanding the non- 
unitarity to second order in small phase shifts A(f> as 
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/(0+A,/)) = /(</))+5.A</)+|A</).iJ.A0+O(A<^), (9) 



with g as the gradient (slope) and H as the Hessian (cur- 
vature) of the non-unitarity. 

The gradient g has 25 components 



df 



(<5no - 5mo) '^{hk - (5j7c) sin(*) (10) 



J>1 



All these Hadamards are drawn from the matrices of 
equi-modular elements {si\ E{(f>) \sm) = e'"^^"' /\/N which — 4>irn + 

we take to be a function on a 25-dimensional phase space, 
with (f) having the individual phases as components: 



2nd column 3rd column 



»26, ■ 



•066), (6) 



6th column 



with sin(*) taking the same argument as before: ★ = 
in — 4'im + 0jm ~ 0jn- It is subjcct to the chain rule for 
differentiation and thus gives rise to the Kronecker deltas. 
The gradient also vanishes for Hadamards because they 
minimize the non-unitarity. 

The 25 X 25 components of the Hessian H are 



9V 



= -E(^^ 



no - Smo){Snp - <5mp) ^((^ifc " Sjk){Sil - Sjl) COs(*). 
j>i 



(11) 



r 



The Hessian is a real symmetic matrix whose eigenval- 
ues give the curvature of /(0) in the directions of their 
corresponding eigenvectors (principal axes). 



Whenever points to a Hadamard, both / and g are 
zero. Moving away by a small step A0 changes the gra- 



dient according to its first order Taylor expansion, 

g{$+A$)^g{$)+H{$)-A$+0{A^). (12) 

To preserve unitarity we require no change to 5 = 
(so that (j) (/) + A<j) continues to minimize the non- 
unitarity). We are therefore primarily interested in the 
nuUspace of H (those vectors n for which H ■ n = 0). 
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IV. NUMERICAL METHODS AND RESULTS 

We calculate numerically^^ the non-unitarity /, gra- 
dient and Hessian H as functions of (f)\ from H we 
calculate its eigenvalues and eigenvectors hi. When 
points to a Hadamard, with the exception of Tao's Ma- 
trix S'g, we always find a 4-dimensional nullspace: four 
eigenvectors hi point in directions of effectively-zero cur- 
vature, i.e. they have numerical eigenvalues Ai < 10""*. 
We keep the four eigenvectors hi that span the nullspace 
and discard the rest. 

From the Fourier matrix F we already knew of four 
directions to move away, namely the (pi corresponding to 
the directions of increasing the parameters (pi of the two 
2-parameter Fourier families, e.g. 

= ( 1,0,1,0,1 , 0, 0, . . . , 0, 1,0,1,0,1 , 0, 0, . . . , 0)/ V6. 
2nd column 5th column 

(13) 

Although the <j)i span the nullspace, they are not all mu- 
tually orthogonal or compatible between the two families; 
moving a finite distance in the (/)2 -I- (/is direction, for ex- 
ample, quickly loses the requisite unitarity; see Figure 

m 

The existence of the nuUspaces, however, suggests that 
a series of infinitesimal moves might preserve unitar- 
ity. The problem is that the nullspace changes from 
Hadamard to Hadamard so that a direction constrained 
to the nullspace at one point might depart from the 
nullspace at another. We propose to "feel a way forward" 
along an evolving direction 9 obtained by projecting the 
previous step's direction into the new nullspace: 

4 

1=1 

with c as a step index. Each step adjusts the relative 
phases ^ by a length in that direction: 

^[c+i] = 0[c] + A6'%]- (15) 

At each step we also use Newton's method to correct for 
any higher order drift (O(A0^)) off the g — Q valley floor: 

4>[c+i\ (f'lc+i] - gic+i]/H[c+i], (16) 

with g and H as the projections of the gradient and Hes- 
sian outside the nullspace; at this point we recalculate 
the Hessian and its nullspace. These corrective steps are 
small (never longer than 1% of AO) and allow us to take 
finite, rather than infinitesimal, small steps = 0.001 
(radians). 

In this way we numerically integrate a parameterized 
curve of Hadamards through the space of matrices of 
equi-niodular elements. Beginning with specifying the 
Fourier matrix F, for example, and taking a small step 
in the (f)2 + (t>3 direction, for example, we improve it with 




FIG. 1: Non-unitarity f23{9) in the cf>2 + (/>3 direction, from 
Fourier; small stepping within the nuUspaces improves uni- 
tarity by 10 orders of magnitude (see scatter plot 10^°/jv(S)). 
Lower panel: the numerically integrated phase shifts. Staying 
within nuUspaces causes a spreading of the phases from their 
original rates of change 0, l/\/T2, and l/\/3. 

Newton's method, obtain a new direction from the pro- 
jection of the previous, take another small step, improve 
it with Newton's method, and so on, progressing a dis- 
tance 9 while maintaining the unitarity, numerically dis- 
covering this curve of Hadamards originating from F and 
the initial choice of direction; see Figure [1] 

As we integrate curves in this manner, we similarly 
evolve four directions 9i spanning the evolving nullspace, 
beginning with, in our example, the four (pi from the 
two Fourier families. In any case, some initial choice 
of nullspace spanning vectors 9^^ can be made. Then 
with each small step along a curve we obtain new direc- 
tions from the previous vectors' projections into the new 
nullspace, 

4 

just as we do for the curve's evolving direction 9[c]- 

The evolving four directions 9i[c] suggest a way to in- 
tegrate Hadamards of order 6 according to four distinct 
parameters 9i: we numerically integrate a distance 9i 
in the first evolving direction, then a distance 02 
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FIG. 2: Four-way foray from Fourier. Upper panel: the nu- 
merically integrated phase shifts obtained by integrating a 
distance 9i — 0.23 in the first evolving direction, &i[c]; then 
02 = 0.17 in the second; Os = 0.34 in the third; and 64 = 0.26 
in the fourth. The first two parameters correspond to moving 
within the first Fourier family; the last two parameters move 
away from it. Lower panel: negligible non-unitarity /]v(S) 
as a function of integrated distance 6 along the four joined 
curves. 



in the second evolving direction, and so on, obtaining a 
Hadamard an integrated distance 6 = J^i along the 
curve of four joined segments. 

In our example, beginning from the Fourier matrix and 
with the 4>i giving the four initial directions, we obtain a 
4-parameter set of numerically integrated Hadamards 

i^6(ei, 02,^3,^4). (18) 

This is our main result. The first two parameters specify 



Hadamards within the first Fourier family. 



(19) 



and the second two parameters then specify Hadamards 
obtained by numerically integrating away from -FgC^i, ^2) 
and away from the first Fourier family in general; see 
Figure [H 

V. SUMMARY: FOUR WAYS FROM FOURIER. 



The 6-level system is the smallest A^-level system for 
which a tomographically complete and efficient set of 
A'' -I- 1 MUBs has not been found; there is a widespread 
search for more than 3. A complete classification of the 
complex Hadamards of order 6, representing bases un- 
biased with respect to a standard basis, would greatly 
narrow the search. The previously known families of 
these Hadamards are topologically connected, perhaps 
by a conjectured 4-parameter family. 

Here we prescribe the numerical integration of a 4- 
parameter set of Hadamards within the space of ma- 
trices of equi-modular elements. We Taylor expand to 
second order a measure of their non-unitarity and use 
the 4-dimensional nuUspace of the curvature to choose 
small steps that preserve the requisite unitarity of the 
Hadamards. We also update, or evolve, a set of four 
directions, as we move along a curve of Hadamards, by 
projecting them into the evolving nuUspace. 

From the Fourier matrix F we choose four initial di- 
rections (two from the first Fourier family and two from 
its transpose) and numerically integrate a 4-parameter 
set of Hadamards Fq{9i, 62, 03,04) by integrating, in or- 
der of increasing z, distances 6i in their corresponding 
evolving directions ^^^c] • 

These results should facilitate numerical searches for 
at least 4 > 3 MUBs in only 12 parameters (4 for each 
candidate MUB beyond the standard basis). They might 
also help point the way to the complete analytical clas- 
sification of Hadamards for N = 6. In that case, then, 
the MUB problem for 6-level systems might finally be 
resolved. 
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